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Understanding  thermal  transport  from  nanoscale  heat  sources  is 
important  for  a  fundamental  description  of  energy  flow  in  materials, 
as  well  as  for  many  technological  applications  including  thermal 
management  in  nanoelectronics  and  optoelectronics,  thermoelectric 
devices,  nanoenhanced  photovoltaics,  and  nanopartide-mediated 
thermal  therapies.  Thermal  transport  at  the  nanoscale  is  fundamen¬ 
tally  different  from  that  at  the  macroscale  and  is  determined  by  the 
distribution  of  carrier  mean  free  paths  and  energy  dispersion  in 
a  material,  the  length  scales  of  the  heat  sources,  and  the  distance 
over  which  heat  is  transported.  Past  work  has  shown  that  Fourier's 
law  for  heat  conduction  dramatically  overpredicts  the  rate  of  heat 
dissipation  from  heat  sources  with  dimensions  smaller  than  the 
mean  free  path  of  the  dominant  heat-carrying  phonons.  In  this  work, 
we  uncover  a  new  regime  of  nanoscale  thermal  transport  that  dom¬ 
inates  when  the  separation  between  nanoscale  heat  sources  is  small 
compared  with  the  dominant  phonon  mean  free  paths.  Surprisingly, 
the  interaction  of  phonons  originating  from  neighboring  heat 
sources  enables  more  efficient  diffusive-like  heat  dissipation,  even 
from  nanoscale  heat  sources  much  smaller  than  the  dominant  pho¬ 
non  mean  free  paths.  This  finding  suggests  that  thermal  manage¬ 
ment  in  nanoscale  systems  including  integrated  circuits  might  not 
be  as  challenging  as  previously  projected.  Finally,  we  demonstrate 
a  unique  capability  to  extract  differential  conductivity  as  a  function 
of  phonon  mean  free  path  in  materials,  allowing  the  first  (to  our 
knowledge)  experimental  validation  of  predictions  from  the  re¬ 
cently  developed  first-principles  calculations. 
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Critical  applications  including  thermoelectrics  for  energy 
harvesting,  nanoparticle-mediated  thermal  therapy,  nano¬ 
enhanced  photovoltaics,  and  thermal  management  in  integrated 
circuits  require  a  comprehensive  understanding  of  energy  flow 
at  the  nanoscale.  Recent  work  has  shown  that  the  rate  of  heat 
dissipation  from  a  heat  source  is  reduced  significantly  below  that 
predicted  by  Fourier’s  law  for  diffusive  heat  transfer  when  the 
characteristic  dimension  of  the  heat  source  is  smaller  than  the 
mean  free  path  (MFP)  of  the  dominant  heat  carriers  (phonons  in 
dielectric  and  semiconductor  materials)  (1-6).  However,  a  com¬ 
plete  fundamental  description  of  nanoscale  thermal  transport  is 
still  elusive,  and  current  theoretical  efforts  are  limited  by  a  lack  of 
experimental  validation. 

Diffusive  heat  transfer  requires  many  collisions  among  heat 
carriers  to  establish  a  local  thermal  equilibrium  and  a  continuous 
temperature  gradient  along  which  energy  dissipates.  However, 
when  the  dimension  of  a  heat  source  is  smaller  than  the  phonon 
MFP,  the  diffusion  equation  is  intrinsically  invalid  because  pho¬ 
nons  move  ballistically  without  collisions.  The  rate  of  nanoscale 
heat  dissipation  is  significantly  lower  than  the  diffusive  prediction 
such  that  smaller  heat  sources  appear  increasingly  inefficient 
in  dissipating  heat.  Furthermore,  heat-carrying  phonons  in  real 


materials  have  a  wide  distribution  of  MFPs,  from  several  nano¬ 
meters  to  hundreds  of  microns  (7).  For  a  given  nanoscale  heat 
source  size,  phonons  with  MFPs  shorter  than  the  hot-spot  di¬ 
mension  remain  fully  diffusive  and  contribute  to  efficient  heat 
dissipation  and  a  high  thermal  conductivity  (or  equivalently,  a 
low  thermal  resistivity).  In  contrast,  phonons  with  long  MFPs 
travel  ballistically  far  from  the  heat  source  before  scattering, 
with  an  effective  thermal  resistivity  far  larger  than  the  diffusive 
prediction.  Phonons  with  intermediate  MFPs  fall  in  between; 
here,  heat  transfer  is  quasiballistic  with  varying  degrees  of  re¬ 
duced  contributions  to  the  conduction  of  heat  away  from  the 
nanoscale  source.  Most  work  to  date  explored  the  reduction  in 
heat  transfer  from  functionally  isolated  microscale  and  nano¬ 
scale  heat  sources  (1-5).  Indeed,  characterizing  heat  transfer 
from  micro/nanostructures  with  varying  size  can  be  used  to 
experimentally  measure  cumulative  phonon  MFP  spectra  of 
materials  (2,  8,  9),  with  the  proof-of-principle  demonstrated  for 
long-MFP  (>l-|im)  phonons  in  silicon  (3). 

In  this  work,  we  show  through  both  experiment  and  theory  that 
the  size  of  the  heat  source  is  not  the  only  important  scale  that 
determines  nanoscale  heat  dissipation.  We  identify  a  new  regime 
of  thermal  transport  that  occurs  when  the  separation  between 
nanoscale  heat  sources  is  smaller  than  the  average  phonon  MFP. 
Surprisingly,  we  find  that  when  phonons  from  neighboring  heat 
sources  interact,  more  of  them  dissipate  heat  in  a  diffusive 


Significance 

A  complete  description  of  nanoscale  thermal  transport  is  a 
fundamental  problem  that  has  defied  understanding  for  many 
decades.  Here,  we  uncover  a  surprising  new  regime  of  nano¬ 
scale  thermal  transport  where,  counterintuitively,  nanoscale 
heat  sources  cool  more  quickly  when  placed  close  together 
than  when  they  are  widely  separated.  This  increased  cooling 
efficiency  is  possible  when  the  separation  between  nanoscale 
heat  sources  is  comparable  to  the  average  mean  free  paths 
of  the  dominant  heat-carrying  phonons.  This  finding  suggests 
new  approaches  for  addressing  the  significant  challenge  of 
thermal  management  in  nanosystems,  with  design  implications 
for  integrated  circuits,  thermoelectric  devices,  nanopartide- 
mediated  thermal  therapies,  and  nanoenhanced  photovoltaics 
for  improving  clean-energy  technologies. 
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regime,  thus  counteracting  the  inefficiency  caused  by  ballistic 
effects  in  the  case  of  isolated  heat  sources.  This  collective  behavior 
can  increase  heat  transfer  to  near  the  diffusive  limit.  Most  im¬ 
portantly,  the  appearance  of  this  “collectively  diffusive”  regime 
mitigates  the  scaling  problems  for  thermal  management  in 
nanoelectronics,  which  may  not  be  as  serious  as  projected  (2, 
10,  11).  Finally,  we  use  this  previously  unobserved  phenomenon 
to  extract  the  contribution  to  thermal  transport  from  specific 
regions  of  the  phonon  MFP  spectrum,  opening  up  a  new  ap¬ 
proach  for  thermal  transport  metrology  and  MFP  spectroscopy. 
This  is  because  by  varying  both  nanostructure  size  and  separa¬ 
tion,  an  effective  phonon  filter  is  introduced  that  suppresses 
specific  sections  of  MFP  contributions  to  thermal  conductivity. 
We  compare  our  extracted  phonon  MFP  spectra  with  pre¬ 
dictions  from  first-principles  calculations  and  find  excellent 
agreement  between  experiment  and  theory.  Looking  forward,  we 
have  a  unique  capability  for  characterizing  phonon  transport  in 
novel  materials  where  predictions  do  not  yet  exist. 

Fig.  1  illustrates  the  differences  between  the  three  regimes  of 
heat  transport  from  nanoscale  heat  sources — purely  diffusive, 
quasiballistic,  and  collectively  diffusive.  Quasiballistic  transport 
(Fig.  IB)  dominates  when  the  size  of  isolated  nanoscale  heat 
sources  is  smaller  than  dominant  phonon  MFPs.  In  the  new 
collectively  diffusive  regime  we  uncovered  (Fig.  1C),  the  sepa¬ 
ration  between  heat  sources  is  small  enough  that  long-MFP 
phonons,  whose  contribution  to  heat  dissipation  would  normally 
be  limited  by  the  small  size  of  nano-heat  sources,  can  once  again 
play  a  significant  role  and  restore  heat  transfer  efficiency  to  near 
the  diffusive  limit.  Although  these  phonons  travel  ballistically 
away  from  each  individual  heat  source,  they  can  scatter  with 


Fig.  1.  Nanoscale  heat  transport  is  determined  by  the  interplay  between 
three  length  scales:  the  size  of  the  heat  sources  L,  the  spacing  of  the  heat 
sources  P,  and  the  MFPs  A,  of  heat-carrying  phonons.  Materials  support  a  broad 
distribution  of  MFPs,  represented  here  by  short  (black)-  and  long  (white)-MFP 
phonons.  (A)  When  all  MFPs  are  smaller  than  L,  heat  dissipation  is  fully  diffu¬ 
sive.  ( B )  As  L  shrinks,  long-MFP  phonons  travel  ballistically,  decreasing  the  rate 
of  heat  dissipation  relative  to  diffusive  predictions.  Short-MFP  phonons  remain 
diffusive.  (C)  When  both  L  and  P  shrink,  long-MFP  phonons  originating  from 
neighboring  heat  sources  interact  as  they  would  if  they  originated  from  a  sin¬ 
gle,  large  heat  source,  enabling  more  efficient  diffusive-like  heat  transfer. 


phonons  originating  from  a  neighboring  heat  source,  thus  creating 
an  effectively  larger  heat  source  size.  In  the  limiting  case,  the 
spacing  between  heat  sources  vanishes  and  this  regime  approaches 
heat  dissipation  from  a  uniformly  heated  layer. 

Experiment 

In  our  experiment,  arrays  of  nickel  nanowires  were  fabricated  by 
e-beam  lithography  and  lift-off  techniques  to  form  periodic  gratings 
on  the  surface  of  sapphire  and  silicon  substrates.  The  nanowire  line 
widths  L  range  from  750  nm  down  to  30  nm,  with  period  P  =  4L 
and  a  rectangular  profile  height  of  —13.5  nm.  The  use  of  nano- 
patterned  structures  rather  than  optical  absorption  allows  us  to 
explore  the  dynamics  of  heat  sources  much  smaller  than  the  dif¬ 
fraction  limit  of  visible  light.  The  metallic  nanowires  are  heated  by 
a  25 -fs  pump  pulse  centered  at  a  wavelength  of  800  nm.  The  sap¬ 
phire  substrate  is  transparent  at  this  wavelength,  whereas  the  silicon 
substrate  has  such  a  long  absorption  depth  that  any  small,  uniform 
heating  of  the  substrate  can  be  neglected.  Laser  excitation  thus 
creates  an  array  of  nanoscale  hot  spots  (lines)  on  the  surface  of  a 
cold  substrate.  Because  all  nanostructures  are  fabricated  on  the 
same  substrate  at  the  same  time,  the  intrinsic  thermal  boundary 
resistivity  at  the  interface  between  the  metallic  line  and  the  sub¬ 
strate  will  be  constant  across  all  samples:  any  variation  in  efficiency 
of  heat  dissipation  as  the  hot-spot  size  or  spacing  is  varied  can  thus 
be  attributed  to  different  regimes  of  thermal  transport. 

The  laser-induced  thermal  expansion  and  subsequent  cooling 
of  the  nanogratings  is  probed  using  coherent  extreme  UV  (EUV) 
light  centered  at  a  wavelength  of  29  nm,  created  by  high  harmonic 
up-conversion  of  an  800-nm  Ti:sapphire  laser  (12).  The  time  delay 
between  the  EUV  probe  pulse  and  the  laser  pump  pulse  is  ad¬ 
justed  using  a  mechanical  delay  stage  between  -400  and  +8,000 
ps,  with  step  size  as  small  as  1  ps.  As  the  EUV  light  diffracts  from 
the  periodic  array  of  Ni  nanowires,  expansion  and  cooling  of  the 
nanogratings  changes  the  diffraction  efficiency,  and  this  signal  is 
recorded  by  a  CCD  camera  as  a  function  of  delay  time  between 
pump  and  probe  pulses.  Examples  of  this  dynamic  signal  are 
shown  in  Fig.  24.  [Note  that  more  details  of  the  experimental 
setup  are  described  elsewhere  (13,  14)  and  in  SI  Text ,  section  SI.] 
Because  the  reflectivities  of  these  materials  do  not  change  with 
temperature  at  EUV  wavelengths  (15),  the  change  in  the  dif¬ 
fraction  signal  can  be  uniquely  attributed  to  physical  deformations 
in  the  surface  profile.  Moreover,  because  the  thermally  expanded 
nanowires  will  change  the  diffraction  signal  as  long  as  their  tem¬ 
perature  is  higher  than  their  initial  state,  this  serves  as  a  measure 
of  the  cooling  rate  of  the  nanowires  rather  than  the  rate  of  heat 
transfer  across  a  specific  distance  as  in  optical  transient  grating  (6) 
or  buried-layer  heating  experiments  (16).  Thus,  the  data  of  Fig.  24 
can  be  used  to  directly  extract  the  average  thermal  expansion  and 
relaxation  of  each  individual  nanowire  induced  by  laser  heating 
and  subsequent  heat  dissipation  into  the  substrate,  in  addition  to 
the  surface  deformations  caused  by  acoustic  waves  launched  by 
the  initial  impulsive  expansion  (13,  14,  17). 

Theory 

To  understand  the  different  regimes  of  thermal  transport  illus¬ 
trated  in  Fig.  1,  we  consider  three  models:  (i)  the  model  described 
in  our  previous  work  that  assumes  isolated  heat  sources  (1);  ( ii )  an 
analytical  model  we  develop  here  to  account  for  interactions  of 
phonons  originating  from  neighboring  heat  sources  using  a  gray, 
single-phonon-MFP  approximation;  and  (iii)  a  more  advanced 
interacting  model  that  includes  a  distribution  of  phonon  MFPs.  As 
discussed  in  detail  below,  this  interacting  multi-MFP  model  allows 
us  to  extract  MFP-dependent  contributions  to  thermal  conductivity 
for  MFPs  as  short  as  14  nm  for  the  first  time  (to  our  knowledge), 
providing  data  that  can  be  directly  compared  with  predictions  from 
first-principles  density  functional  theory  (DFT). 

To  quantify  the  deviations  from  diffusive  heat  transport,  we 
first  build  upon  methods  similar  to  those  described  by  Siemens 
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Fig.  2.  Effective  thermal  boundary  resistivities  are  extracted  from  dynamic  EUV  diffraction.  04)  Dynamic  diffraction  from  60-nm-wide  nickel  lines  on  sapphire 
(Top)  and  silicon  ( Bottom )  display  a  sudden  rise  due  to  impulsive  thermal  expansion  following  laser  heating,  a  long  decay  due  to  thermal  relaxation,  and 
oscillations  due  to  surface  acoustic  waves.  Dashed  black  lines  plot  the  diffusive  prediction,  which  significantly  underestimates  the  thermal  decay  time.  Green 
lines  plot  the  decay  using  a  best  fit  to  the  effective  thermal  boundary  resistivity.  ( B )  Extracted  effective  resistivities  for  each  line  width  L  on  both  substrates 
increase  with  decreasing  line  width  until  the  periods  (equal  to  4 L)  are  comparable  to  the  average  phonon  MFP.  For  smaller  periods  (spacing),  the  effective 
resistivity  decreases  and  approaches  the  diffusive  limit  (black  dashed  line).  The  error  bars  represent  the  SD  among  multiple  datasets  for  the  same  line  width 
samples.  Dotted  red  lines:  predictions  for  isolated  heat  sources  based  on  the  gray  model.  Dash-dot  blue  lines:  gray  model  including  the  onset  of  the  col¬ 
lectively  diffusive  regime.  Solid  purple  lines:  more  complete  model  that  includes  contributions  from  multiple  phonon  MFPs. 


et  al.  (1),  but  including  more  comprehensive  finite-element  physical 
modeling  (18)  to  improve  data  reduction  accuracy.  We  model  our 
system  using  diffusive  heat  conduction  theory,  while  allowing  the 
effective  thermal  boundary  resistivity  (which  sets  the  temperature 
discontinuity  across  the  boundary  between  the  nickel  nanowires  and 
the  substrate)  to  vary  as  a  function  of  line  width  to  account  for 
nondiffusive  effects  in  the  substrate  near  the  heat  source. 

We  use  accurate  sample  dimensions  (height,  line  width,  and 
period)  characterized  by  atomic  force  microscopy  (SI  Text ,  section 
SI).  Fresnel  optical  propagation  is  then  used  to  calculate  the  dif¬ 
fraction  signal  from  the  simulated  surface  deformations.  The  ef¬ 
fective  boundary  resistivity,  reff,  that  provides  the  best  fit  to  the 
experimental  data  of  Fig.  24  represents  the  sum  of  two  terms:  first, 
the  constant  intrinsic  thermal  boundary  resistivity,  rTBR ,  that  origi¬ 
nates  from  the  material  difference  between  nickel  and  substrate; 
and  second,  a  correction  term,  rCorn  due  to  nondiffusive  transport  in 
the  substrate  close  to  the  heat  source  when  either  L  or  P  is  smaller 
than  MFPs.  Although  this  is  a  similar  effective  correction  for  non¬ 
diffusive  transport  near  isolated  nanoscale  heat  sources  as  that  used 
in  previous  works  (2-6),  by  assigning  the  nondiffusive  contribution 
to  the  thermal  boundary  resistivity  rather  than  to  changes  in  the 
thermal  conductivity  of  the  substrate,  we  maintain  a  simple  mod¬ 
eling  geometry  and  avoid  the  need  to  assume  a  particular  region  of 
the  substrate  in  which  an  effective  conductivity  change  should  ap¬ 
ply.  Further  discussion  can  be  found  in  SI  Text ,  section  S2. 

The  effective  resistivity  results  are  plotted  in  Fig.  2 B.  For  large 
line  widths  on  both  sapphire  and  silicon  substrates,  the  effective 
resistivity  converges  toward  a  constant  value — the  intrinsic  ther¬ 
mal  boundary  resistivity.  As  the  line  width  approaches  the  dom¬ 
inant  phonon  MFPs  in  the  substrate,  the  effective  resistivity  rises 
as  thermal  transport  becomes  quasiballistic  and  the  contribution 
to  heat  dissipation  of  long-MFP  phonon  modes  is  suppressed  (1, 
3,  8).  This  behavior  was  successfully  described  in  past  work  using 
a  simple  gray  model  for  sapphire  and  fused  silica,  which  assumes 


a  single-phonon  MFP  to  loosely  describe  a  weighted  average  of 
the  MFPs  from  all  of  the  phonon  modes  contributing  to  thermal 
transport  in  a  given  material.  According  to  this  model,  a  ballistic 
correction  term  proportional  to  Agray/(L/2)  can  be  added  to  the 
intrinsic  thermal  boundary  resistivity  (1,  19);  this  prediction  is 
plotted  in  dotted  red  in  Fig.  2 B. 

However,  as  the  line  width  (and  period)  shrinks  further,  Fig. 
2 B  shows  that  the  effective  resistivity  starts  to  decrease  rather 
than  continuing  to  increase.  The  constant  grating  duty  cycle  for 
our  series  of  samples  means  that  the  smallest-line  width  nano¬ 
wires  are  also  those  with  the  smallest  separation  between  neigh¬ 
boring  heat  sources.  Thus,  for  small  line  widths,  the  separation 
becomes  comparable  to  dominant  phonon  MFPs.  For  silicon,  this 
peak  in  rejf  is  shifted  toward  longer  line  widths/periods  compared 
with  sapphire  because  the  phonon  MFP  distribution  in  silicon  is 
also  shifted  toward  longer  MFPs,  i.e.,  silicon  has  a  longer  average 
MFP  than  sapphire  (1,  6).  As  illustrated  in  Fig.  1C,  in  this  col¬ 
lectively  diffusive  regime,  longer-MFP  phonons  from  neighboring 
heat  sources  interact  with  each  other  as  they  would  if  they  orig¬ 
inated  from  a  single,  large  heat  source,  leading  to  diffusive-like 
heat  dissipation  and  decreasing  the  effective  resistivity.  The 
quasiballistic  model  for  isolated  heat  sources  clearly  fails  to 
capture  this  experimental  observation,  and  a  new  model  for  rCorr 
is  required  to  account  for  the  transition  to  this  previously  un¬ 
predicted  and  unobserved  collectively  diffusive  regime. 

We  propose  to  use  the  concept  of  a  notch  filter  in  the  MFP 
spectrum  to  describe  the  effects  of  grating  line  width  and  sepa¬ 
ration,  shown  schematically  in  Fig.  3/4.  The  notch  filter  sup¬ 
presses  the  contribution  of  phonon  modes  with  MFPs  that  fall 
between  the  line  width  L  and  period  P  of  the  nanogratings.  Thus, 
if  the  grating  period  (separation)  remains  large  while  the  line 
width  is  decreased,  one  would  expect  the  effective  boundary 
resistivity  to  continue  to  rise,  as  shown  in  the  red  dotted  lines  of 
Fig.  2 B.  This  is  because  the  contributions  of  all  phonon  modes 
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Fig.  3.  Line  width  and  period  define  a  suppression  filter  for  phonon  MFP  spectra.  (A)  The  observed  increase  in  effective  thermal  boundary  resistivity  for 
small  line  widths  L  is  due  to  the  suppression  of  the  contribution  to  thermal  conductivity  of  phonon  modes  with  MFP  larger  than  L.  Decreasing  the  period 
P  can  reactivate  modes  with  MFP  larger  than  P,  decreasing  the  effective  resistivity.  In  the  limiting  case  of  a  uniformly  heated  layer,  P  approaches  L  and  all 
phonon  modes  participate  in  thermal  transport.  We  use  as  an  example  the  smoothed  differential  conductivity  distribution  for  silicon  (top  graphs,  green 
line),  calculated  from  first-principles  DFT  (SI  Text,  section  S4).  ( B )  A  comparison  of  the  thermal  decay  of  small  line  width  gratings  for  two  different  periods 
directly  validates  the  prediction  of  the  suppression  filter  model,  i.e.,  small  line  widths  spaced  far  apart  (red  lines)  exhibit  a  slower  initial  thermal  decay 
than  small  line  widths  spaced  closer  together  (blue  lines).  The  dashed  lines  provide  a  guide  to  the  eye  for  the  thermal  decay  through  the  center  of  the 
acoustic  oscillations. 


with  MFPs  longer  than  the  line  width  L  are  suppressed  in  the 
quasiballistic  regime  of  isolated  heat  sources.  On  the  other  hand, 
if  the  grating  period  shrinks,  long-MFP  phonon  modes  start  to 
contribute  again  because  phonons  originating  from  neighboring 
heat  sources  interact  with  each  other  as  they  would  in  a  bulk 
system,  so  the  effective  boundary  resistivity  should  recover  to¬ 
ward  the  bulk  value,  as  seen  experimentally  in  Fig.  2 B.  If  this 
picture  is  true,  then  a  nanograting  with  a  small  line  width  and 
large  period  (i.e.,  wide  filter)  should  exhibit  a  slower  thermal 
decay  than  the  same  line  width  grating  with  shorter  period  (i.e., 
narrower  filter).  Using  a  set  of  periodic  nickel  nanowires  on 
silicon  with  varied  duty  cycle,  we  directly  observed  this  pro¬ 
nounced  difference  as  shown  in  Fig.  3 B. 

To  build  an  analytical  expression  for  rCorr  based  on  this  idea, 
we  use  the  concept  of  a  phonon  conductivity  suppression  func¬ 
tion,  S(L ,  P,  A),  similar  to  those  described  by  others  (8,  20).  This 
suppression  function  is  applied  to  a  bulk  differential  conductivity 


spectrum  versus  phonon  MFP  of  the  substrate,  &(A*),  to  calculate 
an  effective  nanoscale  conductivity  Knano : 

K-nano  =  y.i(Ai)-5(L,P,Ai).  [1] 

We  then  relate  rCorr  to  the  change  in  conductivity  represented  by 
this  suppression: 

rCorr=A{^o~K^k)’  [21 

where  A  collects  geometrical  constants  (discussed  further  in  SI 
Text ,  section  S3)  and  Kbuik  is  the  bulk  conductivity  of  the  sub¬ 
strate,  simply  given  by  A/ ).  For  a  given  phonon  MFP  Ah  S 

must  approach  unity  in  the  diffusive  regime  when  both  L  and  P 
are  large  and  at  the  limit  of  uniform  heating  when  L  =P.  For  the 
limit  of  small,  isolated  heat  sources  when  L  ->  0  but  P  is  large, 
S  L  1(2  A)  to  reproduce  the  behavior  of  the  previously  published 
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quasiballistic  model  (1,  21).  Finally,  the  effects  of  L  and  P  should 
be  the  same  but  opposite  to  each  other  so  that  L  suppresses 
phonon  mode  contributions  in  the  same  way  as  P  reactivates 
them.  To  capture  these  behaviors  along  with  smooth  transitions 
among  regimes,  the  two  effects  are  represented  by  a  special  case 
of  the  generic  family  of  logistic  functions,  where  the  total  sup¬ 
pression  function  is  written  as  follows: 


S  =  tanh 


+ 


1  -  tanh 


[3] 


More  details  of  this  suppression  function  can  be  found  in  SI 
Text ,  section  S3.  Although  more  rigorous  methods  of  deriving 
suppression  functions  for  various  experimental  geometries  are 
currently  being  explored  (9,  20,  22),  none  has  yet  sought  to 
account  for  closely  spaced  heat  sources  or  an  accompanying 
reintroduction  of  phonon  modes.  Eq.  3  represents  the  first  at¬ 
tempt  (to  our  knowledge)  to  include  the  contribution  of  heat 
source  spacing  and  offers  a  model  that  is  simple  enough  for  fast 
integration  into  existing  models  of  heat  transfer  in  nanoscale 
devices,  for  example,  but  complex  enough  to  capture  the  pre¬ 
viously  unobserved  behavior  and  make  successful  predictions. 

To  test  this  new  model  for  rCwr,  we  first  assume  the  simple 
single-MFP  (gray)  model  (where  the  MFP  distribution  is  a  delta 
function).  The  resulting  predictions  are  shown  in  the  blue  curves 
in  Fig.  2 B.  Specifically,  rCorr  in  this  case  is  given  by  the  following: 

rCorr,gray(L,P)  ——  (  ^ 

Kbulk  \S  [T  •>  P •>  Agr(ly  J  J 


Fitting  this  interacting  model  to  the  reff  data  for  sapphire,  we 
extract  values  for  rTBR  and  A gray  that  are  consistent  with  previous 
results  (1):  Agray>  int  =  131  ±  11  nm,  rTBR>  int  =  2.58  ±  0.19  x  1(T9 
m2  KAV.  This  good  agreement  with  the  previous  larger-line 
width  data  and  the  accurate  fit  for  the  full  range  of  our  data 
validate  our  interacting  rCorr  model  as  an  improved  method  to 
account  for  nanoscale  size  effects  in  heat  transport — for  both 
quasiballistic  and  collectively  diffusive  regimes.  Interestingly,  this 
single-MFP  model  provides  a  reasonable  approximation  for  the 
entire  range  of  heat  transport  in  sapphire. 

For  silicon,  the  interacting  rCorr  follows  the  general  shape  of 
the  data  and  yields  Agmy>  int  =  306  ±  17  nm,  which  is  consistent 
with  previously  reported  values  (6).  However,  the  interacting 
gray-model  approach,  although  more  successful  than  the  isolated 
model,  fails  to  capture  the  additional  structured  tail  in  effective 
resistivity  that  appears  for  very  small  line  widths  and  periods, 
below  L  =  100  nm.  The  failure  of  this  approach  is  not  sur¬ 
prising,  because  the  single-MFP  model  is  known  to  be  a  poor 
approximation  for  silicon  with  its  broad  distribution  of  phonon 
MFPs  (8,  20). 

Discussion  of  Collectively  Diffusive  Regime  and  Extension  to 
Phonon  Mode  Conductivity  Spectra 

Having  developed  and  validated  the  new  model  to  capture  the 
transitions  from  diffusive,  to  quasiballistic,  to  collectively  diffu¬ 
sive  regimes,  we  can  now  extend  our  calculations  beyond  the 
simple  single-MFP  model  and  use  our  data  to  extract  the  MFP- 
dependent  contributions  to  thermal  conductivity  in  different 
materials  down  to  MFPs  as  small  as  14  nm  for  the  first  time  (to 
our  knowledge).  Because  line  width  and  period  set  the  location 
and  width  of  the  effective  notch  filter  in  the  phonon  MFP 
spectrum,  each  configuration  uniquely  samples  the  contribution 
to  thermal  conductivity  of  different  MFP  ranges  of  phonon 
modes  with  a  resolution  controlled  primarily  by  the  number 
of  configurations  tested.  The  larger  the  resistivity  correction 
needed  for  a  given  nanograting,  the  stronger  the  conductivity 
contribution  of  phonon  modes  that  were  suppressed. 


Phonon  mean  free  path  (nm) 


Fig.  4.  By  fitting  reff  with  multiple  bins  of  phonon  modes,  the  weights  k( Aj) 
assigned  to  those  bins  give  the  average  relative  contribution  to  the  differ¬ 
ential  thermal  conductivity  (purple  shading).  Both  differential  (distributions) 
and  cumulative  (lines)  conductivities  are  normalized  to  the  total  bulk  con¬ 
ductivity.  For  sapphire  (Top),  our  data  (solid  purple  line)  and  first-principles 
DFT  calculations  (dashed  green  line)  indicate  there  are  no  significant  con¬ 
tributions  from  long-MFP  phonons,  so  the  cumulative  curves  approach  unity 
at  1  |im.  For  silicon  (Bottom),  our  data  are  consistent  with  large  contributions 
from  longer  MFPs. 


To  extract  information  about  the  differential  conductivity 
spectrum,  we  use  the  full  multi-MFP  form  of  rCom  given  by  the 
following: 


rCoJL'P)  A(j2lk(M)S(L,P,  A,)  J'  [5] 

We  partition  the  full  sum  in  bins  according  to  the  MFP  sensitivity 
of  each  grating  configuration  (SI  Text ,  section  S3).  Then  by  fit¬ 
ting  our  set  of  rejf  data  as  rejf  =  rTBR  +  rCorr ,  we  extract  the  average 
k(At)  that  is  associated  with  all  modes  A,  within  each  bin,  thus 
assessing  the  relative  contributions  to  thermal  conductivity  of 
each  region  of  the  phonon  MFP  spectrum  (plotted  in  Fig.  4). 
By  limiting  the  number  of  bins  to  be  no  more  than  one-half  of 
our  number  of  data  points,  we  ensure  a  conservative,  well-con¬ 
ditioned  fit,  although  we  note  that  changing  the  bin  number  does 
not  substantially  alter  the  trends  we  observe.  As  shown  by  the 
purple  curves  in  Fig.  2 B,  for  sapphire,  this  procedure  matches  the 
experimentally  measured  reff  as  well  as  the  gray  model  for  inter¬ 
acting  heat  sources.  For  silicon,  this  more  complete  multiple-MFP 
interacting  model  is  able  to  match  our  experimental  measure¬ 
ments  of  thermal  boundary  resistivity  over  diffusive,  quasiballistic 
and  collectively  diffusive  regimes,  including  the  exceptionally  slow 
drop  in  resistivity  for  small  line  widths  below  100  nm. 

Although  the  number  of  experimental  data  points  limits  the 
number  of  regions  we  can  reasonably  consider  in  this  first  dem¬ 
onstration,  this  approach  still  offers  unprecedented  experimental 
access  to  the  differential  thermal  conductivity  contributions  of  pho¬ 
nons  with  different  MFPs  and  for  benchmarking  theoretical  pre¬ 
dictions,  including  those  from  first-principles  DFT  calculations 
shown  in  Figs.  3 A  and  4.  In  particular,  our  experimental  data  across 
all  MFP  ranges  measured  in  silicon  are  in  very  good  agreement 
with  our  DFT  calculations  [which  also  agree  with  those  in  the 
literature  (8)].  However,  some  small  discrepancies  appear  for 
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phonon  MFPs  around  ~100  nm,  where  experimental  verification 
was  not  possible  before.  Differences  between  the  experimental 
and  theoretical  spectra  in  this  region  may  also  be  exaggerated  by 
our  limited  set  of  small-line  width  gratings;  increased  resolution 
with  a  larger  sample  set  can  address  this  issue.  Our  data  are  also 
consistent  with  observations  by  others  that  emphasize  significant 
contributions  from  long-MFP  (>l-pm)  modes  (2,  3,  6),  but  the 
limited  number  of  data  points  we  have  for  structures  much  larger 
than  the  average  phonon  MFP  results  in  a  relatively  large  un¬ 
certainty  in  this  region.  For  the  purpose  of  comparison  in  Fig.  4, 
we  normalized  the  experimental  spectra  in  silicon  by  assuming 
the  integrated  conductivity  up  to  1  pm  should  match  that  cal¬ 
culated  by  DFT. 

For  sapphire,  both  calculation  and  experimental  data  imply 
that  phonons  with  MFPs  shorter  than  1  pm  are  responsible  for 
>95%  of  the  thermal  conductivity.  The  discrepancy  below  300  nm 
between  experimental  and  theoretical  curves  (apparent  in  the 
cumulative  distribution)  is  due  to  two  factors.  First,  the  sharper  rise 
in  the  DFT  cumulative  curve  arises  due  to  the  very  strong  short- 
MFP  peak  in  the  conductivity  spectrum — a  peak  that  lies  at 
~5  nm,  below  the  lower  bound  of  our  experimental  sensitivity 
(14  nm)  using  30-nm  structures.  Thus,  the  experimental  data  simply 
do  not  include  the  shortest  phonon  MFP  peak.  Second,  because 
of  the  complex  crystal  structure  of  sapphire,  the  DFT  calcu¬ 
lations  required  the  use  of  relatively  small  interaction-distance 
cutoffs  for  determining  the  harmonic  and  anharmonic  force 
constants,  which  may  cause  a  larger  error  in  the  predictions  than 
for  silicon. 

There  are  two  significant  advantages  of  our  EUV-based 
technique  compared  with  previously  reported  MFP  spectros¬ 
copy  techniques.  First,  our  approach  that  combines  nanoheaters 
with  the  phase  sensitivity  of  short-wavelength  probes  is  the  only 
way  to  experimentally  access  dimensions  far  below  100  nm  to 
directly  resolve  the  contributions  of  phonons  with  MFP  down  to 
14  nm.  Other  approaches  require  numerical  extrapolation  tech¬ 
niques  and  interpretation,  which  are  still  being  developed  (9). 
Second,  characterizing  phonon  MFP  spectra  by  harnessing 
interacting  nanoscale  heat  sources  allows  us  to  probe  arbitrary 
segments  of  the  MFP  spectrum  for  any  novel  material,  enabling 
direct  access  to  the  differential,  rather  than  only  the  cumula¬ 
tive,  MFP  distribution. 

The  ability  to  experimentally  extract  a  phonon  MFP  distribu¬ 
tion  down  to  such  small  MFPs  offers  an  exceptional,  useful 


method  for  validating  existing  first-principles  predictions  across 
the  whole  range  of  phonon  MFPs  significant  for  heat  conduction, 
as  well  as  the  first  access  to  such  information  for  more  complex 
materials  where  calculations  have  not  yet  been  performed. 
Furthermore,  combined  knowledge  of  both  the  differential  and 
cumulative  spectra  may  offer  intriguing  insight  into  the  full  MFP 
spectrum  with  the  detail  necessary  for  accurate  prediction  of  heat 
transfer  in  nanostructured  systems. 

Conclusion 

In  summary,  we  uncover  a  new  regime  of  nanoscale  thermal 
transport  that  dominates  when  the  separation  between  nanoscale 
heat  sources  is  small  compared  with  the  dominant  phonon  MFPs. 

We  also  present  a  new  approach  for  characterizing  the  relative 
contributions  of  phonons  with  different  MFPs  that  participate 
strongly  in  heat  conduction.  In  particular,  we  can  probe  the  small- 
MFP  region,  which  has  been  previously  inaccessible  to  experiment. 

This  unique  capability  is  important  as  the  need  for  precise  pho¬ 
non  MFP  distributions  in  complex  nanostructures  becomes  more 
pressing — for  both  fundamental  understanding  and  to  harness 
systems  where  modeling  does  not  yet  exist.  In  the  future,  because 
bright  soft  X-ray  high  harmonic  sources  can  now  reach  wave¬ 
lengths  below  1  nm  (23),  this  approach  can  be  extended  even  fur¬ 
ther  into  the  deep  nanoregime.  In  addition,  more  comprehensive 
and  fundamental  insight  into  nanoscale  thermal  transport  may  be  y 
possible  by  adopting  the  framework  created  to  bridge  all  types  of  g 
anomalous  diffusion  (24-27),  which  does  not  rely  on  effective 
diffusion  models  for  inherently  nondiffusive  heat  transfer.  Fi¬ 
nally,  the  efficient,  collectively  diffusive  regime  of  thermal  trans¬ 
port  that  we  observe  for  the  first  time  can  potentially  mitigate 
projected  problems  for  thermal  management  in  nanoelectronics, 
where  the  power  density  is  likely  to  increase  as  the  individual 
nanostructures  shrink  in  size  (10,  11).  It  also  highlights  important 
design  implications  for  nanostructured  materials  and  devices  for 
energy  and  biomedical  applications. 
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